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Anomalously large thermopower of mesoscopic normal-metal/superconductor interferometers has been in- 
vestigated by Chandrasekhar et al. It was shown that, depending on the geometry of the interferometer, the 
thermopower is either symmetric or antisymmetric periodic function of the magnetic flux. We develop a de- 
tailed theory of the observed thermoelectric phenomena in the framework of the non-equilibrium quasiclassical 
approach. In particular we provide, for the first time, a possible explanation of the symmetric thermopower 
oscillations. This effect is attributed to the electron-hole symmetry violation that originates in the steady-state 
charge imbalance between different arms of the interferometer. Our theory can be tested by an additional con- 
trol over the charge imbalance in a modified setup geometry. We also predict a sign reversal behavior of the 
thermopower with increasing temperature that is consistent with the experiments by Parsons et al. 
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I. INTRODUCTION 

A linear response of a metal to the temperature gradient 
VT and external electric field E is characterized by four 
transport coefficients: a,K,r),C> which are defined, under the 
open circuit conditions, by the following relation 1 

(£Hn)(-vr). <>> 

where J e is the electric current and J® is the heat flow in- 
duced in the circuit. The diagonal terms, a and k, are the 
electric and thermal conductivity, while the off-diagonal ones, 
Tj and (, are the thermoelectric coefficients that fulfill the On- 
sager relation ( = Tr). 

The thermoelectric response is usually observed in a con- 
fined geometry by measuring the voltage rather than the elec- 
tric current. The ratio of the measured voltage to the tem- 
perature gradient applied across the sample is known as the 
Seebeck coefficient (or, the thermopower) and is given by 
S = rj/a = {£ /VT)|je =0 . The Mott's formula^ relates the 
thermopower to the energy dependent conductivity a(e) 

Tr 2 k%T d\na{e) 



and illustrates an important role of the electron-hole asymme- 
try for the thermoelectric response in metals. Indeed, one can 
see from Eq. (0 that S vanishes if a(e) is symmetric around 
the Fermi energy ep. In the case of a spherical Fermi surface 
the thermopower S ~ k^T/eEF is finite only because of a 
tiny difference in the effective masses of electrons and holes. 

In superconductors the situation is complicated by the pres- 
ence of Bose condensate of Cooper pairs that does not respond 
to a temperature gradient. It was suggested by Ginzburg 3 on 
the basis of the two-fluid model of superconductivity that the 
linear response theory (Q3 still applies to a normal (dissipa- 
tive) component of the electron liquid. The thermoelectric ef- 
fect is, however, shunted by the superconducting component 
and cannot be observed in bulk superconductors. 4,5 An incom- 
plete cancellation of the thermoelectric and supercurrents was, 



nevertheless, predicted for anisotropic superconductors and 
superconducting bimetallic samples. 6-7 This theory has been 
experimentally confirmed by Zavaritskii 8 by using a bimetal- 
lic loop to detect the thermopower in a superconducting state. 
The method relies upon the quantization of the superconduct- 
ing condensate in the loop that prevents a complete cancella- 
tion of the thermoelectric current and generates a small, but 
measurable, magnetic flux $t ~ ^0~ 2 ^o, where $o = 2e//i 
is the flux quantum. In subsequent experiments, 9 however, 
much greater values of the thermoelectrically induced flux, 
$ T 10 2 $o 5 were detected. The origin of this "giant flux" 
puzzle remains debated, l&ii 

Mesoscopic systems provide an altogether different way to 
probe the thermoelectric phenomena in the presence of su- 
perconductivity. The superconducting correlations can pen- 
etrate the normal-metal part of the system due to a pro- 
cess known as Andreev reflection. As the result the ther- 
mopower of a mesoscopic normal-metal wire can be affected 
by contacting the wire to a superconductor. First experiments 
of this kind have been performed a decade ago by Chan- 
drasekhar et al. 12 (see Fig. [TJ. In these and in numerous sub- 
sequent experiments 13 ' 14 ' 15,16 ' 17 ' 18,19,20 ' 21 the thermopower of a 
normal-metal wire in proximity to a superconducor/normal- 
metal loop was found to oscillate as a function of the mag- 
netic flux piercing the loop. The amplitude of the oscilla- 
tions was shown to exceed the thermopower of the normal- 
metal wire in the absence of proximity-induced superconduc- 
tivity. The experimental results suggest that the proximity ef- 
fect is responsible for the electron-hole symmetry violation in 
these systems, which is no longer suppressed by a small fac- 
tor fc^T '/ep. The lack of the electron-hole symmetry in the 
experiments was attributed long age 22,23 to a voltage induced 
between the superconducting and the normal-metal part of the 
mesoscopic circuit under non-equilibrium conditions. This 
voltage is directly related to the shift between the chemical 
potential in the normal metal and that in the superconductor. 
Similar shift, which is commonly referred to as the charge im- 
balance, can exist inside a bulk superconductor between the 
chemical potential of quasiparticles and that of the Cooper 
pairs. 

The charge imbalance in bulk superconductors and its 
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role in the thermoelectric effects were intensely studied in 
1970s (see Ref. 24 for the review.). The imbalance can be 
achieved either by injection of a dissipative current to the 
superconductor 2 5,26 ' 27 ' 28 or by applying a temperature gradi- 
ent in the presence of a superconducting flow 29,30,31 (Pethick- 
Smith effect). The relaxation of the charge imbalance is deter- 
mined primarily by the electron-phonon interaction. From the 
detailed theory provided by Schmid and Schon 32,33 one can 
roughly estimate the imbalance relaxation time in bulk s-wave 
superconductors as tq ~ (fcgT/A)Tj n , where A is the value 
of the superconducting order parameter and r in is the temper- 
ature dependent electron-phonon scattering time. In our study 
we deal instead with the voltage drop between the supercon- 
ducting and the normal-metal arms of the interferometer. We, 
however, keep using the term imbalance to stress the analogy 
to the thermoelectric effects in bulk superconductors. 

The non-equilibrium quasiclassical theo ry of p roximity ef- 
fect has been successfully applied in Refs. I34ll35ll36l 37 to de- 
scribe the antisymmetric magnetic field dependence of the 
thermopower in Fig.QJ. In this paper we develop a general an- 
alytical approach that explains both antisymmetric and sym- 
metric dependence of the thermopower observed in the exper- 
iment. Our theory is based on the quasiclassical kinetic equa- 
tions in the limit of weak proximity effect. Even though this 
approximation assumes a small transparency of the normal- 
metal/superconductor interfaces in Fig. [2] it does not affect 
the symmetry of the obtained results. 

In the Section|lV]we solve the kinetic equations for the par- 
allelogram interferometer depicted schematically in Fig. [3] In 
accordance with earlier studies we find that the charging of the 
superconducting arm of the interferometer is determined by 
the interplay of the supercurrent and the temperature gradient. 
The formation of the charge imbalance is analogous to that in 
the Pethick-Smith effect. 29 The extra charge accumulated by 
the superconductor is, therefore, a purely antisymmetric func- 
tion of both the temperature gradient and the applied magnetic 
flux. In this case the left-right symmetry of the structure has to 
be broken in order to observe the effect of proximity-induced 
correlations on the thermopower of the normal-metal wire. 34 
The symmetry breaking can be caused either by a difference 
between L\ and L2 in Fig. [3] assuming that these distances 
are much smaller than the phase-coherence length L^, or by 
a difference in the NS interface transparencies a± and a.2, 
which are given by the ratio of the normal wire resistance per 
unit length to the interface resistance. In the latter case the 
phase coherence of quasiparticles in the normal-metal wire 
is irrelevant for the symmetry breaking. We also stress that 
the energy relaxation processes in the superconductor play an 
important role in the theory of the thermoelectric effect. In- 
deed, the time-independent solution to the kinetic problem ex- 
ists only if the detailed balance condition at a given energy 
is broken. Another words, in the steady-state limit, only the 
total charge transmitted through the superconducting wire is 
conserved while the energy density of the charge flow is not. 
The steady-state regime develops at times exceeding the im- 
balance relaxation time tq inside the superconducting wire. 
The resulting behavior of the thermopower in the parallelo- 
gram interferometer is illustrated in Fig. [4] as a function of 



temperature. 

In the Section [V] we apply the same approach to the house 
interferometer, which is shown schematically in Fig. [5] In 
sharp contrast to the previous case we find a finite contribu- 
tion to the thermopower in the absence of a charge imbalance 
between the normal-metal and the superconducting branches. 
Such contribution has an antisymmetric dependence on the 
magnetic flux and requires a left-right asymmetry of the wire 
N' that must originate in a difference between d% and e?2 rather 
than in a difference between the NS interface transparencies 
Si and a.2. This asymmetric dependence of the thermopower 
in the house interferometer is a phase-coherent phenomena 
that can be attributed to the interference of quasiparticle tra- 
jectories shown in Fig. [6] Apart from the evident requirement 
d\ , da -C L$ this effect is strongly sensitive to the position of 
the Andreev reflection. 38 In the parallelogram interferometer 
the interference contribution to the thermopower is averaged 
to zero since it involves the summation over the trajectories of 
different lengths. 

The odd oscillations of the thermopower are, however, not 
seen in experiments with the house interferometer. Instead, 
a sign-definite thermopower, which has an even dependence 
on the magnetic flux $, has been observed ] 12 ' 16 ' 21 Moreover, 
the maximal values of the symmetric thermopower in Fig.QJ 
correspond to $ = n$o, where n is an integer number and 
$0 = h/2e is the flux quantum. We demonstrate that this be- 
havior is a signature of the finite charge imbalance maintained 
between the Cooper-pair chemical potential in the supercon- 
ductor and the quasiparticle chemical potential in the normal- 
metal. Our theory predicts that the thermopower in this case is 
proportional to fj,(l + cos27r$/$o)i hence its magnetic field 
dependence agrees with the experimental data 12 in FigQJ. The 
thermopower determined by this effect has a peculiar sign- 
reversal behavior with increasing temperature that is illus- 
trated in Fig. [7] This temperature dependence is in qualitative 
agreement with the experiments by Parsons et al. 16 

To our opinion the experiments of Refs. 12.16.21 provide 
a clear signature of a finite steady-state imbalance \i in the 
house interferometer. We demonstrate that a tiny chemical 
potential imbalance between different branches of the house 
interferometer has a great impact on the thermopower due to 
the proximity effect. The possible origin of the tiny imbalance 
is in the usual thermoelectric effect caused by a heat dissipa- 
tion in the second normal-metal wire N'. The validity of this 
scenario can be tested by varying the chemical potential in N' 
with the help of an additional electrode or by a detailed mea- 
surement of the temperature dependence of the thermopower 
in the house interferometer. 

The remainder of the paper is organized as follows. In 
the Section [XT] we introduce the well-known quasiclassical ap- 
proach to the proximity effects in diffusive mesoscopic wires. 
In the Section|III]the rigid boundary conditions at a three ter- 
minal contact are given. The Sections IIVI fVl are devoted to 
the analytical calculation of the thermopower for the paral- 
lelogram and the house interferometer, correspondingly. We 
summarize our conclusions in the SectionlVIl 
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FIG. 1: The experimentally observed thermopower oscillations 
(solid line) in the parallelogram (a) and in the house (b) interferome- 
ter. The figure is adopted from Ref.fT^. 



II. QUASICLASSICAL THEORY 

Andreev reflection is the main microscopic mechanism 
of charge transport between the superconductor (S) and the 
normal-metal (N) at low temperatures. It can be seen from 
the normal-metal side as the conversion of electron-like exci- 
tations of the energy e to the hole-like ones. The reflection 
takes place at the NS boundary and is responsible for the 
phase coherence of electrons and holes at a distance yD/e 
in the diffusive normal metal, where D is the diffusion coef- 
ficient (we set h = kg = 1 throughout the rest of the paper). 
From the superconductor side the same process can be viewed 
as a diffusion of Cooper pairs that brings the superconducting 
correlations into the normal metal. 

The characteristic scale of the proximity effect at a temper- 
ature T is determined by the coherence length £ = y/D/T. 
The electron-hole coherence can be detected in the diffusive 
metal provided T < E c , where E c — D/L 2 is the Thouless 
energy associated with the distance L to the NS interface. On 
the ballistic scales, E c > A, the role of E c is played by A, 
the absolute value of the superconductor energy gap. 

The conditions for the diffusive proximity effect were ful- 
filled in the pioneering experiments with Andreev interferom- 
eters performed by Chandrasekhar et al— and Parsons et ai— 
We, therefore, reduce our consideration to the diffusive ap- 
proximation of the quasiclassical theory that is described by 
the Usadel equation. 



We shall start from the stationary Usadel equation written 
for the quasiclassical Green's function g(x, e) in a diffusive 
normal-metal wire, 



D^ + [iea z ,g] = 0, I 



dx ' 



(3) 



where D is the diffusion coefficient, / is the matrix current, 
and x € {—Li, L 2 ) is the coordinate along the normal-metal 
wire. The Green's function is represented by a matrix in the 
Keldysh space, 



g A 



9 2 = l, 



(4) 



which yields the quasiclassical constraint g 2 = 1, 

Despite the absence of the superconducting order param- 
eter in the metallic wire, the proximity effect will lead to a 
superconducting correlations that are described by the anoma- 
lous components f R ' A of the Green's function. These corre- 
lations are taken into account by the extension to the electron- 
hole space. The spectral (retarded and advanced) sector of the 
Green's function in this space is parameterized as 



hR(A) _(±g R ^ f*^ 

9 ~ y JR{A) Tg R(A) 



(5) 



where bar stands for the charge conjugation. The relation 
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B(A) = J 1 _ fR(A)fR(A) j 



(6) 



follows from the quasiclassical constraint ©. The Keldysh 
component of the constraint, g R g K + g K g A = 0, suggests 
that g K has only two linearly independent entities. In what 
follows we use the standard parameterization 



g K = g R h-hg A , 
where the diagonal matrix 

h = h + a z h„, 



(7) 



(8) 



has two degrees of freedom: h and h a . 

We model the experimental situation by an idealized system 
consisting of two equilibrium reservoirs connected by a non- 
interacting normal-metal diffusive wire N. The distribution 
functions hi and h% in the reservoirs are parameterized by the 
chemical potentials //i, /12 and the temperatures Ti, T2 as 



-Li,iv2 



/ll,2 









where 



h a = tanh ■ 



Ma 



2T„ 



h a = tanh ■ 



Ma 



2T„ 



(9) 



(10) 



and a = 1,2. The functions (1 — h a )/2, (1 — h a )/2 are 
the Fermi distribution functions of electrons and holes, corre- 
spondingly. The electric and the heat current in the normal- 
metal wire are obtained in the Usadel approximation from the 
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Keldysh component of the matrix current /, 

eDu 



J' 



def(s), f(e) = Tra z I K , (11a) 



jQ = x J d££:jQ ( £ ^ = Tri * (llb) 

where I K = g R (dg K /dx) + g K (dg A /dx). 

The Usadel equation (f3]l does not take into account inelas- 
tic processes, hence the charge and heat flow is conserved 
for each energy. Indeed, from the diagonal components of 
its Keldysh sector we obtain the conservation laws 



df(e) _ djQ(e) 



dx 



dx 



= 0, 



(12) 



which play a role of the kinetic equation in the normal-metal 
wire. The diagonal components of the retarded and advanced 
component of the Usadel equation (0 imply the conservation 
of the "spectral" current 



dW i 
dx 



dW / 
dx 



0, 



(13) 



where W R = f R (df R /dx) - f R {df R /dx). This condition 
is equivalent to the conservation of the supercurrent, which 
can flow in the normal metal in the presence of the proximity 
effect. 

From Eqs. (TTTT) we find the current energy densities j e (e), 
j®(s) in the parameterization of Eqs. (I5I7I > as 



dx ax 



(14a) 
(14b) 



where the following functions, which depend on the energy 
and coordinate, are introduced 

M± = (2(1 + .9 V) ± (f R f A + f A f R )) , d5a) 



U = f R f A - f A p 



w = w 1 



w 1 



(15b) 



Thus, the kinetic part of the transport problem is fully de- 
scribed by Eqs. ( 1121141 . 

It is instructive to apply this formalism to the transport in a 
metallic wire in the absence of a proximity effect. In this case 
one finds f R W = 0, hence g R = —g A = a z and 



y dh„ 
dx 



i«=4 



dh 
dx 



(16) 



The solution to Eqs. (fT2l . which yields the boundary condi- 
tions (0, reads 



U \ 1 (i L 2 -L 1 -2x u 
h(e,x) = |(' l oH ^- h T 

h a (e,x) = - ft^ h^ T 



(17a) 
(17b) 



where the following notations are used 

ho = h\ + hi + h 2 + h 2 , 

hx = hi + hi — h-x — h 2 , 

lift = hi-hi+h 2 - h 2 , 

h^r = hi-hi-h 2 + h 2 . 



(18a) 
(18b) 
(18c) 
(18d) 



Substituting Eqs. ( fTTI i to Eqs. ( fT6l ) we obtain j e = 
-(2/L)/i MT andi« = -(2/L)/i T . 

If a temperature gradient is applied between the reservoirs 
that are kept at the same chemical potential \i = /ii = fi 2 , the 
integration in Eq. ( 11 1 at gives J e = 0, hence r/ = 0. The heat 
current in the limit \Ti - T 2 \ <C T is given by Eq. ( Illbl ) as 



J Q - -k{T 2 -Ti), 



tt 2 2DvT 



(19) 



where T — (T 1 +T 2 )/2 is the mean temperature in the system 
and k is the thermal conductivity. 

Similarly, if the chemical potential difference is applied be- 
tween the reservoirs that are kept at the same temperature T, 
the integration in Eqs. (TTTl i gives = and 



r 



M2 - Ml 

(-e)i : 



a = 2e 2 Dv, 



(20) 



where a is the Drude result for the electric conductivity. 

We shall stress that the vanishing of the thermoelectric co- 
efficients 77 and ( in the calculation above is a consequence of 
the linearization of the quasiparticle spectrum near the Fermi 
energy ef- Such a linearization is essential for the quasiclas- 
sical approximation and leads to the exact electron-hole sym- 
metry in metals. The proximity effect can, however, break this 
symmetry through the possible energy dependence of M±, U 
and W in Eq. ( TBI that can give rise to a finite thermoelectric 
response in the system. Such thermoelectric phenomena can 
be analyzed within the quasiclassical approximation. 

Since the superconducting arm of the interferometer used 
in experiments exceeds the electron-phonon scattering length 
it can be described by the equilibrium quasiclassical Green's 
function gs- We assume for simplicity that the supercon- 
ducting energy gap drops abruptly to zero in the normal- 
metal/superconductor interface. Even though this assumption 
is not self-consistent and disregards a small suppression of A 
near the NS boundary, it does not affect our conclusions on 
the symmetry of the thermopower oscillations. 

We, therefore, describe the superconductor by the standard 
bulk expressions for the spectral components of gs that are 
given by 



1 



y/(p T i £ )2 + A2 ^ Ae l * is 



—ie Ae lx 



(21) 



where Aexp(ix) is the superconducting order parameter. 

The external magnetic flux $ piercing the loop leads to a 
gradient of the phase x> hence the phase acquires different 
values xi an d X2 at tne opposite ends of the superconducting 
wire. The order parameter phase difference, 



f = Xi~ X2, 



(22) 



5 



is related to the magnetic flux as ip = 27r$/$o- The Keldysh 
component of g$ is parameterized by 



9§ 



(23) 



where hs = tanh(e/2Ts) is the quasiparticle distribution 
function and Ts is the temperature of the superconductor. 
Here we choose the chemical potential of the Cooper pairs 
as zero and set hs<r = 0. For e < A one finds from Eq. (|2"T1 > 
that g R = gig , hence cjg = 0. In this paper we are concerned 
with ultra low temperatures T <C A and , therefore, disregard 
all quasiparticle effects in the superconductor. 

The effective energy scale for the proximity effect in the 
normal metal is given by the Thouless energy E c = D/L 2 , 
where L is the distance from the NS interface. To describe 
the experimentally relevant situation T ~ E c -C A, we 
restrict our consideration to the excitation energies that are 
much smaller than the absolute value of the superconducting 
gap, therefore 



9§ 



9$ 





,-iX 







(24) 



We will see below that our results for the thermopower are 
indeed defined by small excitation energies e < E c , which 
proves the the consistency of this approximation. 

In addition we assume that the proximity effect is weak due 
to a low transparency of the NS interfaces. This assumption 
allows for the linearization of the Usadel equation (O with 
respect to the anomalous component of the Green's function, 



2 fR 



d 2 f 



dx 



= *7 



2 f R 



z 2 = 



_ 2(0 -ie) 



D 



(25) 



where we replace e by e — iO in order to stress the analytical 
properties of f R . The solution for the advanced Green's func- 
tion is obtained by the substitution e — * —e. The conservation 
of the supercurrent in Eq. (fT3l is equivalent in this approxi- 
mation to the standard property of the Wronskian of a linear 
second-order differential equation. The following symmetry 
relations simplify the subsequent analysis, 

f A (e) = f R (s), /*(*) = /*(-*), f A = (/*)*• (26) 



III. BOUNDARY CONDITIONS 

The boundary condition for a junction shown in Fig.[2]is an 
important ingredient of our description of the proximity effect 
in Andreev interferometers. We assume that a contact wire, 
which is either a superconductor in the parallelogram interfer- 
ometer or a normal-metal wire N' in the house interferometer, 
is described by the Green's function g c . The Green's function 
in the main wire N is denoted by g. 

The transparency of the barrier between the contact and 
the normal-metal wire is parameterized by a coefficient a = 
Tb/(-b- It has a dimension of an inverse length, where Tb is 
the barrier transmission probability per channel and £b is an 
effective barrier length, which is of the order of the mean free 
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FIG. 2: The three-terminal junction formed by a contact attached 
to the main normal-metal wire N. The generalized rigid boundary 
conditions at the junction are given by Eqs. d27|28l > in the matrix 
notation. The Green's function g is continuous along the wire iV but 
there is a jump in the Green's function across the tunnel barrier so 
that g 7^ g at the junction. The kinetic component of Eq. d28| l is 
equivalent to Eqs. J30I3U , 



path in N, The coefficient a can also be regarded as the ratio 
of the normal-wire resistance per unit length to the interface 
resistance and has a dimension of an inverse length. 

The generalized "rigid" boundary conditions 39 apply pro- 
vided Tb <C 1 for each transmission channel at the junction. 
We stress that the tunnel barrier in Fig.|2]is placed between the 
normal-metal wire N and the contact wire. As the result the 
green's function g is continuous along the normal metal wire 
while its derivative has a discontinuity that is determined by 
a difference between g and g c . The continuity of the Green's 
function in the N wire is formally written as 



S[9l 



im (^g{x + 5)~ g(x - 5) J 



lim 

<5 



0. 



(27) 



The discontinuity of the derivative of g(x) is obtained from 
the matrix current conservation in the three terminal junction 



(28) 



Note, that the green's function is discontinuous across the tun- 
nel barrier, i.e. g ^ g c at x = xq. One can easily demonstrate 
that Eq. (|28T i is a discrete analog of the Usadel equation ([3]) 
and yields the same set of conservation laws. 

For instance, the spectral part of Eqs. (1271281 1 leads to 

8 [W] Xq = W = W R - W A , (29a) 
W R = -a (f R f c R - f R f R ) , (29b) 

which illustrates the conservation of the suppercurrent at the 
junction. 

In a full analogy with Eqs. ( 1121141 ) the Keldysh sector of 
Eqs. (127128b is given by 



S\j e (e)} X0 = -Tra z [gc,9}t 
= a(h (r -h aT )M + - a(h-h c )U + -^-^W, 



(30) 



5[j Q (e)h 



T A9c,g]x 



a(h-h c )M- + a{h a -h ca )U + — 



l W, (31) 
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FIG. 3: The parallelogram interferometer. Diffusive normal-metal 
wire N of the length L — Li + L2 connects the temperature reser- 
voirs denoted as Ti and T%. The superconducting wire (dark line) 
is contacting N at two junctions separated by a distance d. The 
TVS' interfaces are characterized by the transparency parameters ai 
and a?2 • Only antisymmetric dependence of the thermopower on the 
magnetic flux is allowed by symmetry. The effect requires £\ 7^ £2 

Or Qi ^ 02. 



where the following notations are introduced 

M± = (g R +g A )(g R +g A ) + l(f R +f A )(f R +f A ) 



+ 2 (f R +f A )if R +f A ), 
u = kf R f A -f R f A -f A f R + f A f R ) 



w 



Mf R f R -f R f R - f A f A + f A f A ) 



(32) 



(33) 
(34) 



The Equations ( 13013 II ) relate the discontinuity in the current 
in the normal-metal wire N to the current flowing into the 
contact wire. 



IV. THERMOPOWER IN THE PARALLELOGRAM 
INTERFEROMETER 

The parallelogram interferometer realized in the experi- 
ment of Ref. 12 is depicted schematically in Fig. [3] The ther- 
mopower is measured between the temperature reservoirs T\ 
and T2 that are connected by the normal-metal wire N (Au) of 
the length L = L\ + L2 ~ 1.7/xra. The wire is contacted to a 
long superconducting arm at the points x = — d\ and x = d 2 
with a distance d = di+d 2 ~ 0.2/im between them. The ther- 
mopower measurement in Fig.QJ is performed at the temper- 
ature T — 350mA". The superconducting coherence length 
at this temperature is estimated as £ = ^/D/T ~ 0.54/ito, 
while the phase coherence length is approximately given by 
~ 3.5/im. Thus, the experimental parameters are such 
that d < £. Therefore, if an external magnetic flux $ is 
piercing the interferometer, the supercurrent J s flows in the 
normal-metal wire between the NS junctions. 

In the subsequent analysis we let [i\ = [12 — [i- and cal- 
culate the quasiparticle current Jo between the reservoirs to 
the first order in the temperature difference T\ — T 2 , thus, ob- 
taining the coefficient 77. The thermopower is, then, calculated 
from the relation S = rj/a. In the limit of weak proximity ef- 
fect the Drude result for a can be used. We note that a similar 



approach has been developed in the series of publications by 
Virtanen and HeikiHa, 34 ! 35 36 

In order to keep our formulas compact we define the dis- 
tances £1 and £ 2 between the NS interfaces and the corre- 
sponding reservoirs as £ a = L a — d a , where a — 1,2. The 
important simplification, which enables us to solve the prob- 
lem analytically, is the linearization of the Usadel equations 
(l25T l for the spectral components of the Green's function g(x). 
This approximation is justified since the wire is in a good con- 
tact to the reservoirs with a vanishing proximity effect, 



— L 1 



= f R{A) 



= 0. 



(35) 



We remind that the chemical potential /i in the reservoirs is 
measured with respect to that of the Cooper pairs in the su- 
perconductor, therefore, the temperature equilibrium is char- 
acterized by /i = 0. 

We start by pointing out that the energy relaxation pro- 
cesses are important for the existence of a steady-state regime 
with a finite /1 in the parallelogram interferometer. Indeed, 
if scattering is purely elastic the charge conservation law in 
the supercondutor must hold for each energy. In this case the 
time-independent solution to the kinetic equation does not ex- 
ist. Following other work a 22 ' 34 we relax the energy-resolved 
condition to the conservation of a total charge in the supercon- 
ductor that is given by 



de [S j 



J 



0, 



(36) 



in the notations of Eq. ((27). The charge conservation con- 
dition in this form can be regarded as the equation on the 
imbalance fi. A finite supercurrent J s in the loop leads 
to a non-zero imbalance, which may act as a source of an 
electron-hole symmetry violation in the system. In the anal- 
ogy with the Pethick-Smith effect, 29 the imbalance defined 
by Eq. (f36l > is proportional to the scalar product of the tem- 
perature gradient and the supercurrent. We will see below 
that \i = A(Tx - T 2 ) sin27T$/$ , where A is a dimen- 
sionless coefficient that is mainly determined by the ratio of 
the Thouless energy E c = D/L 2 to the mean temperature 
T = (Ti + T 2 )/2, such that A is maximal for T - E c . 

From Eq. d36i > we express the current in the superconduct- 
ing wire as 



Js = 



16e 



de 5 



J 



J 



-d-L 



(37) 



where a is the Drude conductivity given by Eq. d20b . Further- 
more, the spacial integration of the electric current along the 
normal-metal wire gives 



T _d je 1 
Jo - Z J s + Z 



dx ^~ 

8e 



def. 



(38) 



Thus, the boundary condition d30T > at the NS interfaces to- 
gether with Eqs. ( II U36l37|38t describe the kinetic part of the 
problem. 
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Let us now perform the calculation of Jo to the second or- 
der in the parameters «i 2- In this approximation the anoma- 
lous components of the Green's function, f R ^ A \ are found 
from the linearized Usadel equation ( f25t with the help of the 
boundary conditions d35l l and 



dx 



-o,\e 



1X1 



-<h 



df 1 



dx 



-a 2 e 



■X2 



(39) 



where the coefficients oti and a 2 introduced in the previous 
Section are determined by the NS interface transparency. The 
conditions d39l follow directly from the retarded component 
of Eq. 



With the help of the following notations 

2 sinh z£ a sinh z(L — £ a ) 



F, 



F\2 



ot\u 2 - 



zL sinh zL 
sinh zl\ sinh z£ 2 
zL sinh zL 



(40a) 
(40b) 



we express the anomalous Green's function near the NS in- 
terfaces as 



e lXl F 1 +e lX2 F 12 , (41a) 
e lX2 F 2 +e lXl Fi2, (41b) 
and find the supercurrent density in the normal-metal wire as 

W=i W * X l ( 42 ) 
\ x<£(-di,d 2 ), 

where Wo 



-4L sin ip Im Fi 2 . From Eq. ( 1 14ab we obtain 

dh n 



= h(-d x )W + U 
U = -h(d 2 )W + A6 



dx 
dh a 
dx 



(43a) 
(43b) 



to the second order in a i >2 . 

The charge conservation condition d36l >. which determines 
the voltage difference /i/e between the superconductor and 
the normal-metal wire, can be written as 



L 



de h T W Q 



jde\5 


dh a 


+ 5 


dha- 


dx 


d 2 


dx 



}• (44) 



The left-hand side of this expression is proportional to the 
temperature gradient due to the distribution function hr- It is 
also antisymmetric with respect to the order parameter phase 
difference ip = xi — X2 — 27r<E» / <I>o due to Wo- The right- 
hand side of Eq. ( |44T > is proportional to the imbalance /i. We, 
therefore, come to the conclusion that the steady-state imbal- 
ance is an odd function of the temperature gradient and the 
applied flux. 

From Eq. ( f3Tb we readily find the supercurrent 



Js 



a 

8e 



de 



h 



Wo 



(' 


dh a 




-5 


dha 


J) 




dx 


d 2 




dx 





(45) 



We also simplify Eq. ( I38l l by calculating the spatial integral 
from the first term in the right-hand side of Eq. ( 114a| ) to the 
leading order in a 



'-Li 

Thus, we obtain 



Jo 



hWdx = ^-[h 



e 2 -h 



h T )W . (46) 




dh n 



dx 



L2 

ILL J 2L L 

-Li -Li 



dh a 

dx 



-di. 



dxM. 



+ f> 



(47) 



where the spatial derivatives of h and h a are found from 
Eqs. (Q/7]). 

The functions and h^x introduced in Eq. ( fT~8T > are sym- 
metric with respect to the energy s and antisymmetric with re- 
spect to the imbalance /1, while ho and hx are antisymmetric 
in e and symmetric in fi. Moreover, the last term in Eq. d47l ) 
can be disregarded in the linear response analysis. Indeed, 
the distribution function h^x is of a second order in the tem- 
perature gradient due to the fact that fi found from Eq. ( |44| ) 
is itself proportional to the temperature gradient. We shall, 
therefore, omit all terms containing h^x in our analysis of the 
thermopower for the parallelogram interferometer. We note, 
however, that the last term in Eq. d47l ) is the only one that 
changes sign with increasing temperature. We, therefore, re- 
gard the sign-reversal behavior of the thermoelectric response 
of the parallelogram interferometer observed by Parsons et 
ali& as the second order effect that is beyond the linear re- 
sponse. Unlike other terms in Eq. d47b the last one is finite for 
a left-right symmetric setup and can easily dominate an ex- 
perimental measurement for any finite temperature gradient. 
We postpone the detailed discussion of this contribution to the 
next Section. 

Omitting the terms containing h^x we find that the discon- 
tinuity of the derivative dh a /dx near the NS interfaces is de- 
termined by the boundary conditions (f30b as 



dh„ 



dx 
dh, 



-di 



dx 



Re(Fi 



Re(F 2 



F 12 coaip), 

F12 COS ip), 



(48a) 



(48b) 



where we take advantage of the condition hsa = 0. 

The spatial integrals in Eq. d47b can be expressed through 
the anomalous Green's functions f R,A at the NS interfaces 
with the help of Eq. (l25T l and the integration by parts. With 
the help of Eqs. (1391411 ) we obtain 



- / dxf R f A = Llm 



F x +F 2 



Fn cos <p), (49) 



that leads to 



L 2 



dx(f R f A -f A f R ) = Q. 



(50) 



-Li 



8 



-i 

"5" 

u 




FIG. 4: (Color online) The function Y(T/E C , d/L, n/2) from 
Eq. J61t versus the ratio T/E c for the parallelogram interferome- 
ter with the symmetric arms £\ = £2- The thermopower is maximal 
for T ~ 2E C , where _B C = D/L 2 is the Thouless energy associ- 
ated with the normal-metal wire N. At high temperatures the ther- 
mopower decays according to the stretched exponential law d60| > that 
is determined by the Thouless energy E' c = D/d 2 associated with 
the distance between the NS junctions. 



Thus, the second term in Eq. d47b is also vanishing. We 
will see in the next Section that this term is responsible for 
the particle-hole interference contribution to the thermopower 
that can be finite in the house interferometer. Finally, the first 
term in Eq. d47b gives rise to the only non-vanishing contribu- 
tion to the current Jo that is linear with respect to the temper- 
ature gradient. 

From Eqs. ( l48b we obtain 



T ad I A h P Fl ~ F2 

Jo = -— / de hp Re . 



(51) 



therefore the quasiclassical thermoelectric effect is absent in 
the left-right symmetric device F\ = F 2 in accordance with 
Ref. [22I and general symmetry considerations. The ther- 
mopower is readily found from the relation S = ( Jo /a) (T 2 — 
Tj.) -1 , where we disregard the corrections to the Drude con- 
ductivity a arising from the proximity effect. This approxi- 
mation is applicable to the second order in the NS interface 
transparency. 

In order to calculate Jo in Eq. (fSTb we substitute // found 
from Eq. ( l44l into Eq. ( 15 It . Disregarding the terms that are 
proportional to h^T we obtain the equation on the imbalance 
/i in the following form 



sin tp — I de fix Im F12 



= j de h„ Re + F 2 ) + F 12 cos ^ . (52) 

Expanding this equation to the first order in /i and the temper- 
ature gradient we obtain 



[i = A(Ti - T 2 ) siny>, 



where T = (Tf + T 2 )/2 is the mean temperature in the system 
and the coefficient A is given by 



.4 



(d/L) J dee cosh ^ ImF 12 
T J de cosfT 2 ^ Re(F x + F 2 + 2F 12 cos <p) ' 



(54) 



We note that the coefficient A is finite even for vanishing in- 
terface transparency parameters 0(1,2- From Eqs. ( 151153b we 
find the thermopower 



Ad sin tp 
8eT 



J de cosh" 2 ^ Re(F 2 - Fx), (55) 



that behaves roughly as sin ip. The thermopower is finite only 
if there is a difference between F\ and F 2 , hence the device 
asymmetry is required. The Equation d55l ) describes the ana- 
log of the Pethick-Smith effect in a superconducting proxim- 
ity system. 

The imbalance given by Eq. d53l is a non-monotonous func- 
tion of the mean temperature T that reaches its maximal value 
at T ~ E c = D /L 2 and decays as a stretched exponent for 
T ^> E' c = D/d 2 . It is worth noting that the imbalance itself 
is finite even for a left-right symmetric setup and in the limit 
of small transparency of the NS interfaces. The temperature 
dependence of fi is determined by that of the coefficient A in 
Eq. ( f54b . The energy integrals in Eq. d54"i i can be rewritten 
as sums over the Matsubara frequencies. This greatly simpli- 
fies the calculation of the asymptotic behavior of A at large 
temperatures. 

Let us estimate the coefficient A for a symmetric setup with 
l\ = l 2 — I and a\ — a 2 . In the low temperature limit 
T < E c we obtain 



.4 



T Li 



1 - 2{e/Lf 



18l-(2£/J) 2 sin%/2) 



(56) 



For high temperatures T ^> E' c we find from the Matsubara 
representation of Eq. ( f54l that 



.4 



where 



Cl 



(2-2- 1 /2)^(3/ 2 ) 



0.93. 



(57) 



(58) 



The analysis shows that the coefficient A reaches its maximal 
value at T ~ E c = D/L 2 . 

For T 3> E' c the energy integral in Eq. d55l ) is decaying 
as T -1 / 2 , provided ot\ 7^ a 2 . This decay crosses over to a 
stretched exponential decay for a\ = a 2 if the asymmetry is 
caused merely by a difference between £1 and £ 2 . Therefore 
an asymmetry in the transparency parameter is of a greater 
importance to the thermopower than that in the arm lengths. 
For an equal arm interferometer £\ = £ 2 = £ with ot\ 7^ a 2 
we obtain from Eqs. (156157b in the limit T < E c = D/L 2 
that 



(53) 



_ dc v {a 2 - a{) soup 



TJ 2 

15 



(59) 
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For high temperatures, T ^> E' c — D / d 2 , we get 



-U 



S 



Trd(a 2 — a\) sin<^ d 2 
8^ 1? 



Ik. 

2ttT 



- v /2nT/E' c 



(60) 

The stretched exponential decay of the Pethick-Smith ther- 
mopower at high temperatures is entirely due to the behavior 
of the imbalance \i. 

For the interferometer with l\ = £2 — £ it is convenient 
to introduce the dimensionless function Y by the following 
relation 



5' 



ird(a 2 — a 2 ) sin (f> 



Y 



T 



d 



(61) 



We see from Eqs. (178160b that Y is a very smooth func- 
tion of the order parameter phase difference <p. The depen- 
dence on <p disappears in the limit T ^> E' c . We plot the 
temperature dependence of Y in Fig. [4] for ip = tt/2 and 
d/L = 1/4, 1/3, 1/2, 2/3, 3/4. The thermopower reaches its 
maximal value in the interferometer with l\= £2 = d = L/3 
at the temperature T ~ 2E C . Our results are consistent with 
the experiment of Ref. and with the numerical analysis of 
Refs. 34.35.40. The magnetic flux dependence of the ther- 
mopower is always antisymmetric and is essentially given by 
sin27r$/$o- A small deviation from this law may arise at 
low temperatures T < E' c . The strength of the proximity in- 
duced thermoelectric effect in the parallelogram interferom- 
eter is restricted by the asymmetry parameter a\ — a 2 and 
by the maximal value of the dimensionless function Y that is 
approximated by Y max w 0.05. 



V. THERMOELECTRIC EFFECT IN THE HOUSE 
INTERFEROMETER 

In this Section we apply the same theory to the house in- 
terferometer depicted in Fig. [5] In our description of the sys- 
tem we introduce two normal-metal wires, N and AT', that are 
connected via a single tunnel junction with the transparency 
parameter a. In the actual experiment of Ref. 12 the junction 
NN' is ballistic because the metallic wires are produced in a 
single lithography circle. We, however, expect that the sym- 
metry of the thermoelectric coefficient is not affected by the 
detailed characteristics of the NN' junction. 

Different space variables x and y are introduced as the coor- 
dinates along the wires N and N', correspondingly. The prox- 
imity effect is described by the Green's functions g c (y) in N' 
and g{x) in N. The N N' junction is located at x = y = and 
the wire N' is contacted to the superconductor at y = —d\ and 
y = c?2. The transparency parameters ai and 0,2 of the N'S 
interfaces are introduced in the same way as in the previous 
Section. 

The thermopower measurements of Ref. [l2] are shown in 
Fig. QJ>. In the experiment the distance between the N'S is 
given by d = di + d*z ~ 2pm that is smaller than the phase 
coherence length L$ ~ i.Qfim at T = 38mK. We, therefore, 
expect that the setup is adequately described by the Usadel 
equation. 




T 2 



T',V 



FIG. 5: The house interferometer. The wires N and N' are normal 
metal wires. The dark wire is superconducting. The transparency 
parameters a, ai,2 are defined by the ratio between the transmis- 
sion probability per channel and the effective barrier length. Both 
symmetric and antisymmetric dependence of the thermopower on the 
magnetic flux $ is possible in this setup. The antisymmetric effect, 
however, requires di 7^ di. A charge imbalance between the super- 
conducting wire and the temperature reservoirs is necessary for the 
symmetric thermopower. 



The charge conservation condition in the loop formed by 
the superconductor and the N' wire is expressed as 



J 



0. 



(62) 



Thus, there exists no Pethick-Smith like contribution to Jo 
that is analogous to the first term in Eq. ( [47) . Another words 
the mechanism of the superconductor charging due to a super- 
current flow is absent in the house interferometer. Unlike in 
the parallelogram interferometer, the other two contributions 
to Jo in Eq. d47l i survive, 



Jn = 



8e 



de < 



h T 1 
2LL 



dxY - 



1 



2L L 



dx M 4 



-hi 



-h. 



(63) 

It is worth noting that the expression d63l for the dissipative 
current between the temperature reservoirs is valid for any in- 
terface transparency. 

The first term of Eq. i63\ , which is proportional to Kt, is 
antisymmetric in the magnetic field and is further referred to 
as the interference contribution. This contribution is insensi- 
tive to the imbalance and vanishes in the parallelogram inter- 
ferometer because of the property ( |50b . It remains, however, 
finite in the house interferometer where the length of the inter- 
fering trajectories is fixed by the distances d\ and c?2 between 
the N'S interfaces and the AW junction. 

Two interfering quasiparticle trajectories, which contribute 
to the heat transfer between the reservoirs, are depicted 
schematically in Fig. [6] The scattering amplitude of the pro- 
cess can be written as 

A e - a 2 \a ie t{k '- kh)dl+lXl + a 2 e l{k '- kh)d2+lX2 1 2 , (64) 

where fc^ = kpzLe/vp are the wave vectors of electrons and 
holes. Here we took into account that the Andreev reflection 
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FIG. 6: An example of interfering trajectories that contribute to the 
thermoelectric response of the house interferometer. 



process is of a second order in the N'S interface transparency. 
The corresponding contribution to the thermopower is propor- 
tional to the energy derivative of the scattering amplitude A £ 
at the Fermi energy, 



S 



de 



ex a 2 a\a 2 {d 2 — d\) s'vcup, 



(65) 



which is a manifestly antisymmetric function of the magnetic 
flux. The result of Eq. d65l l is finite if d\ ^ d 2 , while the 
difference in the transmission parameters a\ and a 2 is neither 
necessary nor sufficient for the observation of the interference 
effect. Integration over all possible trajectories in the paral- 
lelogram interferometer lead to the complete suppression of 
interference that is reflected in Eq. (l50l . 

The interference contribution to the thermopower discussed 
above has not been observed in experiment, because it is over- 
whelmed by the second term in Eq. ( l63b which we refer to as 
the symmetric contribution. This term is proportional to h^T 
and symmetric with respect to the magnetic flux provided \l is 
flux independent. It remains to be finite for a left-right sym- 
metric device but vanishes for fi = 0. 

The experimental measurements of Ref. 12 shown in 
Fig- 01 can be regarded as demonstration of a finite imbalance 
state formed in the house interferometer. The experiment sug- 
gests that a constant imbalance /i is largely independent on the 
temperature gradient in the wire TV as well as on the supercur- 
rent in the loop. Since /x = is the only solution to the charge 
conservation condition d62l ), the origin of such a steady state 
is not entirely clear. We stress, however, that a tiny charge 
imbalance can originate in the usual thermoelectric effect pro- 
vided N' wire is in an equilibrium state with a temperature T' 
that is smaller than both T\ and T 2 due to the phonon cooling. 
In this case we can estimate fi — e(T — T')Sn>, where Sn' is 
the small thermopower of the N' wire. The tiny charge imbal- 
ance /i will be seen to have a great effect on the thermopower 
due the proximity effect enhancement factor. In what follows 
we simply keep /i as a phenomenological parameter of the 
quasiclassical theory. 

We start by calculating the anomalous component f R (x) 
of the Green's function in the normal metal wire N. We again 
take advantage of the linearized Usadel equation (l25T l and em- 
ploy the rigid boundary condition 



df^ 
dx 



(66) 



where = f R (y = 0) is the value of the anomalous Green's 
function at the junction. 

In addition to Eq. d66| ) we use Eq. ( f35T > and the continuity 
of f R (x) at x = 0. As the result we obtain 



f R = 



«/5 



sinh zL 2 sinh z(L\ + x), —L\ < x < 



z sinh zL \ sinhzii sinhz(L 2 -a;), < x < L 2 

(67) 

We calculate the spatial integrals in Eq. d63l with the help of 
Eq. d67| i as 



1 



- / dx (f R f A ± ff) = £ (/M ± f2>f£) Im F, 



f A 7R 



D 



2e 



R 7 A 
cOJcO 



f A 7R\ 



F = a- 



. sinh zLi sinh zL 2 
zL sinh zL 



L = Li+ L 2 . 



(68a) 



It is worth noting that the first term of Eq. (1631 1 is not neces- 
sarily vanishing unlike in the parallelogram interferometer. 

We evaluate the second term of Eq. (l63l l to the leading order 
in the N'S interface transparency by calculating the following 
integral 



1 
L 

K 



dx(f R f R +m 



A IA\ 



D 

2s 



lmf R f R (K-F), 



2 L\ sinh 2 zL 2 + L 2 sinh 2 zL\ 



L sinh zL 



(69) 



that is also expressed through the the anomalous Green's func- 
tion in the N' wire. Thus, from Eq. ( l63l l we find, to the second 
order in / c q, that 



Jo = 



f^f R )ImF 



(70) 



2L 



{(f?of£> + tfofJo) Im F + Im f R f R (F - K )) } . 



It remains to determine the function f c (y) from the Usadel 
equation in the wire N'. Since we calculate the current Jo to 
the leading order in the interface transparency parameters the 
tunneling between and N' can be disregarded in the Usadel 
equation. Taking advantage of the rigid boundary conditions 
in the form 



dfl 
dy 

we obtain 



f R (y)=^^ 



dfl 
dy 



= & 2 e lx \ 



(71) 



<h 



coshz(d 2 -y) . 

a\e Al — — h a 2 e 

z smnza 



lX2 coshz(di + y) 



: sinh zd 



(72) 

where d = d\ + d 2 . Thus, the value of the anomalous Green's 
function at the contact is given by 



c0 



_ t cosh zd 2 _ . coshzdi 
a i e — -,+a 2 e lX2 ^- — -. 
z sinn zd z sinh zd 



(73) 



Substitution of Eq. d73l to Eq. ( |70l completes the calculation 
of the thermoelectric effect. 
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Let us now analyze two different terms in Eq. d70i i sepa- 
rately. We first let /i = so that h^T = and the only contri- 
bution to Eq. (TTOb is antisymmetric in the magnetic flux. From 
Eq. (l73l we find the expression 



f, 



R fA 
cOJcO 



fcofco = 4aia 2 sin</? 



Im (cosh zdx cosh z*d^) 



\z sinh zd\ 2 

(74) 

In order to obtain the thermopower S = (Jq/(j)(T2 — Ti) -1 
we substitute Eq. d74T i to Eq. (T70b and expand to the linear 
order in the temperature gradient as 



hr = — 



sin 



T 2 cosh 2 fi, 



(75) 



The resulting expression for S reads 

a 2 5t\ OLi sin tp Dd 2 B 
~ 8eL f 

with the coefficient B given by 



B 



Im (cosh zd\ cosh z*d'i) tanh(zL/2) 
de n Im ■ 



| zd sinh zd | 2 T cosh 4* zL sinh zL 



(77) 



where we let Li — L 2 = L/2 for simplicity. The Equa- 
tion d76l ) describes the interference contribution to the ther- 
mopower. The result of Eq. ( l76b is indeed proportional to 
a 2 aio;2 sin tp and is vanishing for d\ = d 2 in accordance with 
the estimate fl65l . Moreover, in the limit T <C min {.E c , £^,}, 
where £ c = D/L 2 and = D/d 2 , we find 



5 : 



a a\5i2{d2 — d\)DL 
198eTd 



■ sin tp, 



(78) 



that is equivalent to Eq. II 65b. We note that the divergence 
of Eq. (fTST l in the limit T — > is regularized by the high 
order processes in the parameter a. The result of Eq. d78l is, 
therefore, valid only for T > a 2 D. More accurate expression 
can be obtained by solving the non-linear Usadel equation in 
the normal-metal wire N' . 

For higher temperatures T ^> max {E C ,E' C } it is legiti- 
mate to substitute cosh(e/2T) — > 1 in Eq. ( 1771 ). The remain- 
ing integral gives the dependence of the thermopower on the 
distances L, d\, and g?2- Thus, unlike the Pethick-Smith con- 
tribution d6Tb to the thermopower, the interference contribu- 
tion has a monotonous temperature dependence and decays as 
the power law T~ 2 at high temperatures. It is also the only 
contribution that remains finite for /i = 0. 

Let us consider now the symmetric contribution to Eq. ( T70b 
that emerges for a finite fi. Since this contribution does not 
rely on the device asymmetry, we let for simplicity d\ = c?2 = 
d/2. Following the discussion above we treat the imbalance 
\i as a constant in our quasiclassical analysis and calculate the 
thermopower with the assumption that fi is independent on the 
temperature gradient T% —T-y. The substitution of Eq. d73l l to 
Eq. ( f70b gives 



Jn = 



a D(a\ 



a 2 + 2a.\Ui cos^j) 



4e 



TL 



I 



deh^P, (79) 
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FIG. 7: The function Y(T/E C ) from Eqs. J87l88t versus the ratio 
T/E c for the house interferometer with di = d.2- The function is 
maximal for T ~ 4E C and changes sign at T « 24E C , where E c = 
(76) d/ L 2 

is the Thouless energy associated with the normal-metal wire 
N. The inset shows the temperature dependence of the thermopower 
in Eq. d87t under the assumption that the imbalance fi and the length 
L are temperature independent. The thermopower is maximal at T ~ 



Ec and decays as T~ £ for T > 2AE C 



where 



T 

P = Im 

16s 



F -K 



|zsinhzd/2| 2 2(zsinhzd/2) 2 



(80) 



Thus, for the completely symmetric case ai = &2 the ther- 
mopower has an overall phase dependent factor 1 + cos tp, 
that qualitatively agrees with the experimental curve shown in 
Fig. HJi. In general, the minimal absolute value of the sym- 
metric contribution to the thermopower is determined by the 
asymmetry parameter (<Si — d^) 2 - 

We note, however, that the integral in Eq. ( |80l diverges at 
small energies. This divergence is regularized by the high or- 
der processes in the parameter a. In order to get a well-defined 
expression we need to solve the non-linear Usadel equation 
in N' that is a daunting task. We give instead an analytical 
estimate to the symmetric thermopower by replacing the N' 
wire with a diffusive quantum dot. In this model the Green's 
function, g c , in the dot can be obtained in the spirit of the 
Nazarov's circuit theory 41 from the matrix current conserva- 
tion condition 



[«i5s,i + otig§,2 + aa Zl g^} = 0. 
Resolving the condition for we obtain 

/ c iI = (aie iXl +a 2 e iX2 )u- 1 , 

with 

u = ((Si + a 2 )(-ie/A) + ay/l + (-ie/A) 2 )' 
+ a\ + & 2 + 2ot\a.2 cos (p. 



(81) 



(82) 



(83) 



In this approximation the anomalous Green's function — 
ftfj acquires no spacial dependence on y. The function 
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remains to be small in the parameters a,\ and there- 
fore Eq. ( TTOb is still a legitimate approximation for the ther- 
mopower to the leading order in <Si.2- Substituting f^ Q = 
from Eq. (l82l to Eq. (T70l > we arrive at the result d79l with 



F -K 



T f F 
P = - Im < r-r , 



(84) 



In order to be consistent we have to take the limit e <C A and 
a a *C a in the expression ( f83t , hence 

u = ev 2 , (85) 

and the parameter a cancels out completely in Eq. (1841) . In the 
limit \i < T and T x - T 2 < T we expand /i^t as 



T 2 cosh 2 e/2T 



0-f 



tanh 



27V 



(86) 



and substitute this expression into Eq. d79l ) with the function 
P given by Eq. d84l >. This provides us with the final result 
for the thermopower S = (Jq/<t)(T2 — Pi) -1 of the house 
interferomter 

P>/i(a 2 + u\ + 2&i5 2 cos <p) - / T \ 

where F c = D/L 2 and the dimensionless function Y(T/E C ) 
is given by 



y = 



2a; tanh x — 1 f 3 tanh « 3 
eta; ^ Im • 



4a; cosh x 



1 



1k x 



cosh ft x 



, (88) 

yJ-ixT/E c . In the limit T < E c we find 
(187tt 2 /7560)(P/F c ) 3 . The function Y deter- 



With /Cj; 

Y(T/E C 

mines the temperature dependence of the symmetric contri- 
bution to the thermopower (see Fig. |7]i provided [i is a tem- 
perature independent constant. This contribution is not sup- 
pressed by any asymmetry factor unlike the Pethick-Smith 
thermopower (l6Tb and the interference contribution dToT ). 
Moreover, the symmetric contribution is neither monotonous 
nor sign-definite function of temperature. At T ~ F c the 
thermopower approaches its maximum value, which can be 
estimated as 



5'| 



T~E a 



n 1 A^Va 



2ciici2 cos <p). 



(89) 



Thus, the presence of a very small imbalance n leads to large 
observable magnetic flux dependence of the thermopower. 
The absolute value of the thermopower at T ~ E c and $ = 
is determined by the parameter /j,L/(£bE c ), where £b is the 
effective barrier length that is orders of magnitude smaller 
than L. This is the proximity effect enhancement factor that 
gives rise to an anomalously strong sensitivity of the ther- 
mopower to the charge imbalance /i. 

We also see from Eqs. J87|88t that the thermopower 
changes sign at T w 24P C and decays as T~ 2 at higher 
temperatures since Y — > —2/5 for T ^> 24P C . For higher 



temperatures the length L has to be substituted by the phase 
coherence length provided the latter is smaller than the 
distance between the reservoirs. The temperature dependence 
of Lcj, complicates the direct comparison to the experimental 
data. In particular this dependence should strongly enhance 
the reversed thermopower for T > 24P C . Nevertheless, the 
sign reversal behavior of the the symmetric thermopower ob- 
served by Parsons et ali&H is in a qualitative agreement with 
the result of Eq. d87l 

A thorough experimental test of the presented theory can 
be performed with the help of an independent experimental 
control over the imbalance fi that has a strong effect on the 
sign and the magnitude of the symmetric contribution to the 
thermopower in the house interferometer. For the parallelo- 
gram interferometer the symmetric contribution is described 
by the last term in Eq. j47l i, that is of the second order in the 
temperature gradient. This term can be roughly estimated by 
substituting /i from Eq. (T53l > to Eq. (f8Tb . This second order 
effect reverses sign with increasing temperature and is anti- 
symmetric in $ due to the magnetic flux dependence of the 
imbalance in the parallelogram interferometer. It is likely that 
the effect of this type has been observed in the experiments by 
Parsons et al^^- 



VI. CONCLUSION 

In this paper we applied quasiclassical theory to study the 
origin of the thermoelectric effects in Andreev interferome- 
ters. The theory predicts three additive contributions to the 
electric current expressed by Eqs. dl la|14at in the presence of 
the superconducting proximity effect. Depending on the ge- 
ometry of Andreev interferometer and its parameters any of 
these contributions may dominate the quasiclassical thermo- 
electric response of the normal-metal wire. One can classify 
the observed thermoelectric effect by its dependence on the 
temperature and the magnetic flux piercing the interferome- 
ter. The last term in Eq. (I14ab . which describes the proxim- 
ity induced supercurrent, and the second term, which is re- 
lated to the interference contribution, are antisymmetric in the 
magnetic flux. Even though the supercurrent cannot flow be- 
tween the normal metal reservoirs the last term in Eq. ( 1 14ab 
can, nevertheless, contribute to the thermopower by means of 
the proximity-induced Pethick-Smith effect. This effect takes 
place provided the temperature gradient is aligned in a part of 
the normal metal wire with the supercurrent. This situation is 
realized in the parallelogram interferometer considered in the 
Section HV] 

The left-right asymmetry of the parallelogram interferom- 
eter is responsible for a difference in the dissipative charg- 
ing currents flowing from the normal-metal reservoirs to the 
superconductor. As the result of this asymmetry a com- 
pensating dissipative current Jo dBTI i is flowing between the 
normal-metal reservoirs. The analytical expression for the 
thermopower in the case of weak proximity effect is given 
by Eqs. (|54|55l l. The proximity-induced Pethick-Smith ef- 
fect demonstrates a non-monotonous temperature dependence 
with a maximum at T ~ 2P C (see Fig.|4| and is characterized 
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by the stretched exponential decay for high temperatures ( f60b - 
The effect does not rely on the quasiparticle phase coherence 
in the normal-metal wire. The most important ingredient of 
the theory is the difference in the electron distribution function 
near the first and the second NS interface in the parallelogram 
interferometer. 

In contrast, the interference contribution to the ther- 
mopower, which is given by the second term in Eq. (114al i. 
does require a phase coherence in the normal-metal wire. This 
contribution is the only one that exists in the absence of a 
charge imbalance between the superconductor and the normal 
metal. This effect cannot be seen in the parallelogram inter- 
ferometer but might be observed in the house one provided 
d\ 7^ c?2 in Fig. [5] The interference contribution is character- 
ized by a monotonous temperature decay, that is estimated in 
Eqs. ( |76l77t in the Section|V] 

Finally, the first term in Eq. (I Hal l is regarded as the sym- 
metric contribution. It is largely insensitive to the left-right 
symmetry of the device and gives rise to the thermopower 
that is symmetric with respect to the magnetic flux (with the 
assumption that [i is flux independent). In the Section [Vl we 



argue that the existing experiments indicate the presence of 
such a steady state imbalance in the house interferometer. The 
symmetric contribution is characterized by a peculiar temper- 
ature dependence in Eqs. (187188b that is neither monotonous 
nor sign-definite. This contribution strongly affects the ther- 
moelectric response of the parallelogram interferometer in the 
second order with respect to the temperature gradient. 
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